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Abstract. Dynamical random walk of classical particle in thermodynamically equi- 
librium fluctuating medium, - Gaussian random potential field, - is considered in the 
framework of explicit stochastic representation of deterministic interactions. We dis- 
cuss corresponding formally exact Langevin equations for the particle's trajectory and 
show that Marcovian kinetic equation approximation to them is inadequate, - even 
(and especially) in case of spatially-temporally short-correlated field, - since ignores 
such actual effects of exponential instability of the trajectory (in respect to small 
perturbations) as scaleless low-frequency diffusivity/mobility fluctuations (and other 
excess degrees of randomness) reflected by third-, fourth- and higher-order long-range 
irreducible statistical correlations. We try to catch the latter, - squeezing through 
typical theoretical narrow bottleneck, - with the help of an exact relationship be- 
tween the instability and diffusivity statistical characteristics, along with standard 
analytical d approximations. The result is quasi-static diffusivity fluctuations which 
generally are comparable with mean value of diffusivity and disappear in the limit of 
infinitely large medium's correlation length or infinitely small correlation time only, 
in agreement with the previously suggested theorem on fundamental 1 /f noise. 
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1. Introduction 

In this paper we shall touch classical analogue of interesting quantum statistical- 
mechanical problem already touched in [H |2]. Namely, statistics of dynamical "ran- 
dom walk", or "Brownian motion", of (microscopic) particle interacting with thermo- 
dynamically equilibrium scalar boson field. For instance, with (harmonic) phonon field 
(crystal lattice or other medium oscillations). 



Key words and phrases, dynamical foundations of kinetics, molecular Brownian motion, random 
walks and mobility 1/f fluctuations in (infinitely) many-particle Hamiltonian systems, fundamental 
1/f- noise . 
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In [U|2], basing on the widely known Hamiltonian model of such interaction ( "polaron 
model"), we obtained an exact system of shortened evolution equations for probabili- 
tity distribution of the "Brownian particle" (BP) and its statistical correlations with 
the "phonon field" ("boson thermostat"), and argued that exact solutions to these 
equations includes 1/f-type (scaleless) low-frequency fluctuations of BP's diffusivity 
and thus mobility. But direct formulations of these solutions, or at least good enough 
approximation to them (without loss of the diffusivity/mobility 1/f -noise), still are 
absent. Therefore, for further progress in statistical-mechanical theory of fundamental 
1/f-noise, it may be useful to consider classical variant of the mentioned Hamiltonian 
model. 

At that, we can avoid any detailing of "phonon" (thermostat) part of full Hamilton- 
ian, if instead exploit the "stochastic representation" (SR) of dynamical (deterministic) 
interactions, for the first time suggested and tested in [3j 0] and later generalized, de- 
veloped and applied in [5J El El E] . Thus our consideration will be at once additional 
probing of old and search of new SR possibilities. 



2. Principles of the stochastic representation 



Firstly we have to recall main SR statements [31 HI [7]. Let a Hamiltonian sys- 
tem consists of a "Dynamical subsystem" under our interest, "D", and some its envi- 
ronment, or "thermal Bath" (thermostat), "B", and full Hamiltonian of "D+B" has 
bilinear form 

H = H d + H b + H int , H int = At B n , (1) 

n 

where operators (or phase functions, in classical mechanics) H d , D n and H\, , B n are 
defined in Hilbert spaces (or phase spaces) of "D" and "B", respectively. Besides, let 
initially, somewhen in the past, full density matrix (probability distribution function) 
of "D+B" p{t) had factored form: p( m ^ = pit — > — oo) = p d n ^ p£ . Then marginal 
density matrix (DM) of "D", pd(t) = Tr& p(t) , can be expressed as average value, 

Pd (t) = (P(t)) , (2) 

of randomly varying DM p(t) satisfying stochastic von Neumann (or Liouville) evolu- 
tion equation 

dp(t) 



dt 



h 



pit), H d + J2 X n(t) D Ti 



^2 Vn{t) D n O pit) 



(3) 



where o denotes symmetrized (Jordan) product, A o B = (AB + BA)/2 , and x n (t) 
and y n {t) are random processes. 

At that, all statistical characteristics of x n (t) and y n (t) are unambiguously de- 
termined by internal dynamical properties of "B" , along with its initial DM, . 
Corresponfing formulae can be found in [3j HJ |5j EJ [7] (for most general variants of SR, 
including non-Hamiltonian dynamics, see [5]). In particular, if "B" is a set (contin- 
uum) of harmonic oscillators (wave modes), while p[ m ^ has canonical Gibbs form, with 
some temperature T , then x n (t) and y n (t) are stationary Gaussian random processes 
representing thermodynamically equilibrium "Gaussian thermostat" . 

Below, we confine ourselves by this case, but generalize it to continuous index n : 



J2 D nB n ^ J D{r)B{r)dr , (4) 



^x n {t)D n ^ J x(t,r)D(r)dr , y n (t) D n => J 



y(t, r) D(r) dr 



where r marks points of a d -dimensional space, and dr = d d r . 

3. Random fields of equilibrium thermostat 

It is important to remind, firstly, that x n (t) and hence x(t,r) represent direct 
dynamical perturbation of "D" by "B", i.e. thermostat noise, while y n (t) and y(t,r) 
inverse perturbation of "B" by "D" and related feedback action of "B" onto "D", in 
particular, "friction" or "viscosity", etc., i.e. thermostat induced dissipation. Therefore 
y n (t) or y(t,r) are peculiar ("ghost") random variables: all their self-correlation are 
zeros, e.g. (y(ti, ri) yfe, ^)) = 0, althouh their cross- correlations with x n {t) or 
x(t, r) can differ from zero. 

Secondly, the generalized fluctuation-dissipation relations (FDR , see e.g. [TUl [HI 
[121 [TBI [141 fT5l fl6l ITTl fT8l [19] and references therein), resulting from fundamental 
properties of Hamiltonian dynamics, imply definite mutual correspondence of the "x-y" 
cross-correlators and "x-x" self-correlators . In particular, for equilibrium and spatially 
homogeneous thermostat, according to standard recipes from [31 HIES [6], we can write 

duj 



K xx (r,r 1 -r 2 ) = {x(t + T,r 1 )x(t,r 2 )} = / cos ut S(u 1 r 1 - r 2 ) n 

J-oo 271 

K xy (T,r 1 -r 2 ) = {x{t + T,ri)y{t,r 2 )) = (5) 

2 r . , fa N duo 

sin cur tanh — &{u,ri — r 2 ) 



CO 
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where 9{r) is Heaviside step function, thus expressing both "x-x" and "x-y" correlators 
through one and the same spectral function, S'(c < ;,r 1 ,r 2 ) > 0. Importantly, in accor- 
dance with the causality principle, y(t,r) is only correlated with later x(t' > t,r') . In 
the classical limit, the FDR (El) reduces to 



K xy (r,r) 



e{T) lKUr,r) 



(6) 



T Or 

In case of harmonic thermostat, these two correlators completely determine (Gaussian) 
statistics of x(t, r) and y(t, r) . 

4. Particle in thermal random field 

To consider "Brownian particle" (BP) in a thermally fluctuating media, let us model 
the latter with potential field and choose in (jTJ)- (J4J) 

H d ^- , D n D(r) = 5(r - R) , (7) 
2m 

with P and R being (operators of) BP's momentum and coordinate and B(r) (oper- 
ator of) the fluctuating potential. Then, in the Wigner representation, the "stochastic 
von Neumann (quantum Liouville) equation" ([3]) takes form 



dp P „ _ i 

-£ = — v rP + - 

at m n 



r [t, R - y V P 



ih . 



.r[t,R+j Vp 



ih 



y[t,R-—Vp)+y[t,R+—V 1 



P + 



P 



with \^ and Vp denoting derivatives (gradients). In the classical limit, clearly, it 
turns to "stochastic Liouville equation" (SLE) 



-£ = Vpp + Vx{t,R)V P p + y{t,R)p, 

at m 

where Vx(t, R) = V R x(t, R) . 

5. Langevin equations 

Natural solution to the SLE (j5J) is 

p{t, R, P) = 5(R - R{t)) 5(P - P{t)) exp [ / y(t', R(t')) dt' 

Jt>t> 

where R(t) and P(t) are random processes obeing stochastic ODE 
dR(t) P(t) dP(t) 



(9) 



dt 



rn 



dt 



-Vx(t,R(t)) 



(10) 



1 Notice that in [31 H] the minus sign in §5§ was absent, but this misprint had not penetrated to 
results of SR applications. 



with some (may be random) initial conditions in the past. 

These equations do not display dissipative feedback action of the media (i.e. BP's 
self-action through media), which however is completely involved by field y(t,r) and 
ecomes apparent after averaging expression ( flUl) to get (J2]). 

At the same time, classical SR allows [7] to perform in ( flUj) separate averaging 
over "x-y" cross-correlations, in such way removing y(t,r) from f JTOj) and dispaying 
its effect in ffTTj) . As the result, stochastic equations ( TTTj) transform to what can be 
called "Langevin equations" (LE). A general recipe to construct LE is accumulated 
by Eqs. 39-42 in [7]. In case of Gaussian thermostat, it strongly simplifies (see Eqs.44 
and following example in [7]), and in application to our present system, as defined by 
EqsJH yields 

dRjt) P{t) 
dt m ' 

f(t,R(t)) - f K(t-t',R(t) - R(t'))dt' (12) 
Jt>t> 



dP(t) 
dt 



't>v 

Here, the integral represents medium's feedback response,, introducing friction and 

dissipation, its (vector-valued) kernel K(t, r) is expressed by 

Iff dtu d d k 

K(t,r) = V r K xy (t,r) = — J j sin kr sin ur kuj S(u, \k\) — -^-^ , (13) 

and f(t,r)) = — Vx(t, r) is feedback-free "seed" medium's random force field possess- 
ing Gaussian statistics with zero average value, 

(f(t,r)) = 0, (14) 

and correlation function 

( f(h, ri) * f(t 2 , r 2 ) ) = V ri * V r2 K xx (ti ~ t 2 , n - r 2 ) (15) 

( * denotes tensor product of vectors). Now the random distribution function (DF) 
to be inserted to © is, instead of (P2J, merely pit) = 5(R - R(t)) 5(P - P(t)) . 

Of course, BP's state {R(t), P(t)} is densely correlated with the force f(t,R(t)), 
therefore Gaussianity of the field f(t, r) in itself does not mean Gaussianity of f(t, R(t)) 

6. Quasi-quantum formulation 

Excluding from theory, - in the spirit of quantum mechanics, - BP's momentum, and 
considering BP's coordinate marginal distribution only, 

W{t, R) = [ p{t, R, P) d d P , 



one can derive for it, instead of ©, equations 

dW ^ R) = -V[F(t, R) W{t, R)) + y(t, R) W(t, R) , (16) 

V(t,r) = - —VA(t,r) , 
m 

where V(t,r) is random velocity field generated by scalar "action" field A(t,r) which 
satisfies nonlinear stochastic Hamilton- Jacobi equation: 

^|^ = ^-(VA(t,r)) 2 + x(t,r) (17) 

Then, again performing partial averaging in respect to "x-y" cross- correlations, we can 
transform this into 

dW ~ 1 

= -V[VW] , V = VA , (18) 

at m 

dA (VA) 2 
~dt ~ 2m 



t>v 



+ x(t, r) + / K xy (t -t',r- r') W(t', r') d d r' dt' (19) 



Obviously, the latter stochastic PDE envelopes EqsJT2]and serves as a kind of LE. 

7. CONCRETIZATION OF PROBLEM 

It is seems reasonable to assume, first, that our fluctuating media is statistically 
isotropic, therefore functions from ([5J) are spherically symmetric functions of coordinate 
differences, and we can write 

d d k 



S(u, r) = J cos kr S(u, \k\) j—^ 

(with, clearly, S(—cj, k) = S(lo,k)). 

Second, the force field is not too singular, so that BP's momentum diffusivity is 
finite: 

POO 

/ (f(t,R(t))*f(t-T,R(t-T)))dT ~ 

poo 

~ / [-V* VK xx (r,R(t) -R(t-r))]dr = D {p) ^ oo (20) 
Jo 

This condition holds only if 

k 2 S{Vk,\k\)d d k < oo (21) 



at any velocity V . Then BP's momentum and velocity, V(t) = dR(t)/dt = P(t)/m. 
behave as continuous random processes. 



Third, at any k , 

J S(w, \k\)duo < oo (22) 

This means that f(t,r) possesses nonzero correlation time, r c , i.e. is not a "white 
noise" in time. Evidently, otherwise space-time variations of f(t,r) would be able to 
propagate with infinite velocity, which would be rather nonphysical behavior. 
Fourth, the spectral function S(u>, \ k\) is such that 

lim VK xx (t,Vt)) = 

T->0 

Then, using FDR (JS} and calculating the feedback integral in EqJT2] by parts, we can 
transform EqsJT2l into 

dR(t) 



dt 

dV(t) 



V(t), 

f(t, R(t))/m - I G(t- tf, R(t) - R(t>)) V{t') dt' , (23) 
Jt>t' 

where matrix response function 

G(r,r) = (/(r,r)*/(0,0))/mT = -V *VK xx (r,r)/mT = 

I f f duj 

cos kr cos cut k * k S(lo, \k\) — dk (24) 



mT J J 2tt 
determines BP'viscous friction (here and below dk = d d k/(27c) d ). 

It is interesting task to reveal, under enumerated conditions, long-time statistical 
properties of BP's displacement, e.g. AR(t) = R(t) — R(0) . 

8. Likely conjectures and Marcovian approximation 

8.1. Time-local friction approximation. Notice that because of condition (T221 the 
friction always is more (may be highly) or less (may be slightly) nonlinear in respect 
to BP's velocity. 

If characteristic velocity's relaxation time, 

r v ~ mT/D {p) 

(with Z)( p ) 's diagonal in mind), is large enough in comparison with that of G(r, Vt) 
(i.e. the force's correlation time r c ), then one can visualize the non-linearity by "time- 
local approximation" of LE (|23|) . 

dR(t)/dt = V(t) , 
dV(t)/dt = f{t, R(t))/m - g(V{t)) V(t) , (25) 



with time-local relaxation rate 



g(V) = / G(r,Vr))dr 



2mT 



k*kS(Vk, \k\)dk 



(26) 



8.2. Marcovian approximation. Simultaneously, the above approximation pushes 
to treat f(t,R(t)) like delta-correlated ("white") noise. In order to reasonably de- 
termine its characteristics, we have to consider evolution equation for the random DF 
p(t) = p(t, R, V) , corresponding to EqsfJSl i- e - SLE 



dp/dt = {-VW R + Vv[g{V)V - f(t,R)]}p 



(27) 



and coarsen it into approximate kinetic (Fokker-Planck) equation for actual BP's DF 
pdit) = pdit, R,V) d2D- Quite standard manipulations yield 



d Pd (t)/dt 



-VV R + V v g(V)V] Pd (t) - V v m- l (f{t,R)p{t)) 

-mr x {f(t,R)p{t)) * 

- / (f{t, R) e-^v fit - r, R) V v pit ~ r) ) dr * 
^ Jo 



m 

oo 



(fit, R) e~ TVWR fit-r,R)} V v e rVVR dr 



o 

oo 



(f{t, R) fit - t ,R-Vt)) (Vv + rV R ) dr 



(f(t,R)f(t-r,R-Vr))dr 



Pdit) = 
Pdit) R 
VvPdit) 



that is finally, in the simplest "one-loop" approximation, 

dp d /dt = -VV R p d + V v giV) [iT/m)V v + V\ 



(28) 
(29) 
(30) 
(31) 

(32) 



Thus, we in fact replaced the force f(t,R(t)) by BP's coordinate-independent but 
instead velocity-dependent Gaussian white noise f(t,V(t)) with correlator 

(f(t,V(t))*f(t',V(t'))} v = 2Tmg(V)5(t-t') ~ 2D^5(t-t') , (33) 

with (. . . )y meaning conditional averaging under fixed V. 

8.3. Marcovian stochastic equations. To write out an equivalent SE, notice, in 
view of the g(V) 's definition ( 12 6 p and condition (121 p . that one always can make single- 
valued smooth change of variables, V =>■ U , such that 



(dU/dV) g(V) idU/dVy 



9 , 



9 

where ~g is a constant (unit matrix). It will be good choice if it equals to g(V) 's 
average over equilibrium Maxwell probability distribution of velocity: 

^J 9 (V)M ( V )d *V, MiV) = (34) 

(naturally, - as Eqj32] says, - stationary V 's distribution is Maxwellian) § . Then, in 
terms of U , the SE corresponding to Eqs l32ll3"3l looks merely as 

dR/dt = V(U) , dU/dt = f(t)/m - j(U) , (35) 

7 (t7) = gVu[V 2 (U) - (T/m) In det g{V{U))}/2 , 

with velocity-independent white noise source: 

(f(t)*f(t'))v = 2TmgS(t-t') (36) 

9. Inadequacy of Marcovian approximation: 
disclosing of conventional conjectures 

For the first look, we just demonstrated that our problem hides nothing novel, since 
reduces to quite trivial SE. But this is wrong impression. 

The matter is that the Marcovian approximation ( I32l -(l36"l) has qualitative defect: it 
neglects the above mentioned non-Gaussianity of the force acting onto BP, f(t, R(t)) , 
and therefore losses specific non-Gaussian (higher-order) correlations between f(t, R(t)) 
and BP's path R(t) . This loss took beginning in transition from expression (|28|) to 
expression ( 1291) . which just means replacement of f(t,R(t)) by Gaussian white noise. 
In fact, the resulting Eqs J32ti36l may arouse suspicions already because have no essen- 
tial difference from equations for particle under short-correlated in time but infinitely 
far-correlated (constant) in space random force! 

In order to better feel importance of the loss, let us compare "fidelities" of solutions 
to SE (|35|) and LE ( 125]) . that is their sensibilities to small perturbations (e.g. that 
of initial conditions). At that, non-linearity of friction plays no essential role, and for 
simplicity and visuality we deal with linear friction. 

9.1. Fidelity of solutions to SE in Marcovian approximation. First, consider 
differential response of EqsJ35ls solutions to infinitesimally small change of (initial) 
velocity at time t — , that is 

v(t) = dV(t)/dV(0) , r(t) = dR(t)/dV(0) 

2 Notice also that singularity of such change of variables would indicate inapplicability of time-local 
approximations at all, not speaking about Marcovuan approximation. 
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Since the noise source in Eqs]35]is insensible to BP's state, it disappears under differen- 
tiation in respect to V(0) or U(0) , that is does not influence on v(t) , r(t) . For linear 
friction, when U — V and j(U) = gV , with g = g = const , we thus have equations 



with initial conditions v (0) = 1 , r(0) = . Consequently, at large t velocity's pertur- 
bation certainly tend to zero, while path's (coordinate's) one to a constant: 



9.2. Fidelity of solutions to LE in time-local friction approximation. Now, 
turn to the approximate but more adequate LE (1251) . From them we have, also at 
linear friction, 



again wiyh initial condotions v(0) = 1 , r(0) = . Thus, now we meet essentially multi- 
plicative "noise source", (V/(t, R{t))/m) r , which can not be made state-independent 
by a non-singular change of variabes. On average, solution to these equations coincides 
with (1371) . But in the sense of fluctuations it is much more interesting: naturally, it is 
statistically unstable. 

To see this in most simple way, let us consider ID case, d = 1 , - when v (t) , r(t) , 
etc., become scalars instead matrices. Introduce random DF 



dr/dt = v , dv/dt = — g v , 



t;(*) -)• , r(t) -> 1/g 



(37) 



dr/dt = v 



dv/dt = (Vf(t,R(t))/m)r - gv , 



(38) 



g(t, r, v) = S(r — r(t)) S(v — v(t)) 



and derive approximate kinetic equation for its average 



g(t,r,v) = (g(t,r,v)) 



by treating f(t,R(t)) as white noise (like in (J29l)-(l311). The result is 



dg/dt = [-vV r + gV v v)g + Q{rV v f g , 



(39) 



where parameters are expressed by 




(40) 
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Next, considering evolution of second-order statistical moments, from Eqj39] (or di- 
rectly from Eqs]38]we have 

(r 2 > \ / (r 2 ) \ / 2 \ 

= M (rv) , M = -g 1 (41) 
V (v 2 ) J \2Q -2g J 

Eigenvalues, \ij , of matrix M. are roots of cubic equation 

ti(ji + g)(ti + 2g) = 4Q (42) 

It clealy shows that one of roots, - let be denoted by /x+ , - is real positive, that is 
solutions of Eqsj38] are unstable in the sense of second-order (and hence higher-orfer) 
moments. 

Rate of the instability, compared with velocity relaxation rate g ~ 1 /r v , is 

determined by dimensionless paramer Q/g 3 . From (gQJ) it follows that 

Q T A 2 

g3 mr^g 2 r 2 ' 

where r c is characteristic correlation length of the force field f(t,r), and A ~ 
y^T/m/g ~ r v y^T/m characteristic BP's "free path" length. Hence, if r c C A, 
then Q/g 3 ^ 1 and, according to Eq|42j n + ~^> g . 

9.3. Fundamental incompleteness of Marcovian approach and typical "bot- 
tleneck" of dynamical theory of noises. We just revealed crucial defect of Marco- 
vian approximation: it completely losses exponential instability of BP's trajectories in 
respect to their small perturbations. As the consequence, it losses all statistical effects 
of this instability and therefore, generally speaking, can give only a caricature of real 
dynamical noise. 

Then, how one should avoid the loss? The answer was prompted long ago by critical 
analysis of "molecular chaos" in fluids |20j, crystals [22], under charge transport [191 
|23~1 [24] and generally in transport phenomena J2TJ [23]. Namely, we have to reject 
any a priori statements (even very attractive) about "independencies" of "random" 
variables, - like e.g. "molecular chaos hypothesis" or "marcovianity" , - and allow 
any statistical dependencies and correlations compatible with equations of statistical 
mechanics. After that we might find that some of a priori unexpected or neglected 
dependencies and correlations really take place and are physically important. 

Thus, in general there are two variants of approximate theory of transport noise: 
one (conventional) before the mentioned theoretical "bottleneck" (overcoming usual 
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instinctive conjectures§) and another behind the "bottleneck". Figuratively speaking, 
that are two different solutions of same approximate equations, "trivial" and "non- 
trivial". The latter contains low-frequency fluctuations, - like 1/f- noise, - of transport 
characteristics which are constants in the former^. 

10. First steps to solution of the problem. Exponential instability, 
higher- rder statistics, and scaleless diffusivity fluctuations 

To understand possible consequences of the exponential instability, first it is useful 
to point out several facts. 

10.1. Some useful formulae and remarks. 1. Stationary (equilibrium) distribu- 
tion of BP's velocity has Maxwellian form ( 134"1) . i.e. is Gaussian, regardless of degree of 
BP-thermostat interaction. This statement follows from the structure of Hamiltonian 
of our system as defined by (JT|) and ([7]). Therefore, in equilibrium statistical ensemble, 
for any function $(V, . . . ) of V and some other random factors independent on V , 
we can write 

( V, <3>(V, . . . ) ) =(T/m)(V y $(V,...)> , (43) 
{V,V, $(V,...)) = (T/m) 2 (V v V v $(K,...)> , 

etc. Here and below, angle brackets with n commas denote joint cumulant of n + 1 ex- 
pressions separated by the commas (the Malakhov's cumulant brackets [38]). Similarly, 
since field f(t,r) is Gaussian, 

(f(t,r),<t>(VJ,...)) = Jdt' J dr'(f(t,r),f(t',r')) ( ^^° ) ( 44 ) 

(so-called Furutsu-Novikov formula). 

2. BP's state {R(to), V(to)} at arbitrary chosen (and then fixed) "initial" time mo- 
ment to can be considered as statistically independent on the simultaneous medium's 
state. Then, any function of later BP's states {R(t),V(t)} (t > to) and BP's dis- 
placement (path) AR = R(t) — R(to) gives an example of the mentioned function 

3 "Prejudices" disclosed by N. Krylov in [25] , 

4 For more explanations see review-discussion pages in [TJ [SJ [TU [501 US 1221 ?,IH] and also ?, 
HH [2U [301 ED 1321 [331 1331 [33- A hidden role of the exponential instability in fundamental 1/f- noise 
formation was directly demonstrated [22 . 
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$(V, . . . ) , with V = V(t ) , and Eqf43] yields, in particular, 

( V(t), V(t ) ) = (T/m) ( dV(t)/dV(t ) > , 

( AR, V > = (T/m) ( dAR/dV ) , (45) 
( AR 2 , V , V ) = (T/m) 2 (d 2 AR 2 /dV 2 ) , (46) 

where for brevity we use Vq = V(to) . 

3. Notice that at t — t ^> r v second of expressions f|4"5l) gives BP's diffusivity, let 
be denoted by D . According to EqsfSZl D ~ T/mg = T/rrig. Thus, mean value of 
the differential response dAR/dV is directly connected to the diffusivity. 

Similarly, fluctuations of this response are closely connected to BP's diffusivity fluctu- 
ations. The latter, on the other hand, can be adequately characterized by fourth-order 
BP's path- velocity cumulants. Namely [21 UHl EUJ EH E31 EH EE! EHl E21 EH EH], at 
r = t — t ^> t v , function 

K D (r) = ^^(AR,AR,AR,AR) = 

= lil( A i?( 4 )) = lil(Ai?( 4 )) = -^--f-(AR^) (47) 
24 dt 2 v 1 24 dt§ X 1 24 dtdt^ 

plays role of correlation function of equilibrium BP's diffusivity fluctuations (and that 
of BP's mobility fluctuations in weakly non-equilibrium regime under external force). 
Here ( ) is short designation for X 's n -order cumulant, and we took into account 
that at t — t > t v the path statistics must depend on time difference t — to only. We 
thus have 

K D (t -t ) = ~( V(t) , AR, AR , V(t ) ) , (48) 

or equivalent ly 

K D (t - to) = I ( AR, AR, V(t ), V(t ) )~l(AR, AR, AR, dV(t )/dt ) = (49) 
2 o 

= I ( V(t), V(t) , AR, AR ) + I ( dV{t)/dt , AR, AR, AR ) 
2 6 

The differential response appears in visual form if we notice that at t — 1 ^> t v the 

two velocity values, V{t) and V{t ) , certainly are almost statistically independent 

and hence mutually Gaussian random quantities. Consequently, formula ( 148]) can be 

transformed as follows, 



2 



T 

K D (t-t ) => 



1 / d 2 AR 2 \ _ / OAR \ / OAR 

2 \dV(t)dV(t ) / ~ \ dV(t)/ \ dVito) 



(50) 
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Clearly, right-hand side here consists of squared first-order differential response and 
besides second-order (double- differential) one. Analogously, with the help of Eqs 141111441 
and other above formulas, transforms expressions fj49l) . 

4. Integrating EqJ4Tl one in the same fashion obtains relations 



The latter thus calls for analysis of third-order irreducible correlations (cumulants) 
which address Kd(t to mutual correlation between "fidelity" and "diffusivity" ofBP's 
trajectories. In essence, this is particular case of general exact expression for correlation 
functions of low-signal excess noise and dissipation fluctuations [22] (see also [TO]). 

10.2. How the instability might work. 1. In the Marcovian approximation, the 
fourth cumulant (AR^) certainly is either linear function of t — t — to at r > r v or 
(in linear approximation) identical zero. Hence, Kd(t) can appear non-zero at r > r„ 
due only to what is lost under Marcovian approximation, i.e. the exponential instability 
of BP's trajectories. Hence, to calculate some of the K d {t) 's expressions fj47|) - fl5T|) . one 
have to substitute there exact solutions of LE (fT2"|) or f[2"3"j) and then perform necessary 
averaging over realizations of the field f(t, r) . However, it is impossible in literal sense, 
just because of the exponential instability. 

At present, we do not know a regular method to break away from this "vicious 
circle" (penetrate through the "bottleneck"). Therefore, it would be quite good if we 
demonstrated significance of (AR^) and Kd{t) at t > r„ at least under some 
reasonable approximation of LE. 

In this respect, the last of equivalent expressions fj47|) - fl5T]) appears most suitable, 
in combination with the approximate time-local Eqsf55] and EB1 since this combination 
most visually highlights statistical interference between BP's trajectory {R(t),V(t)} 
itself and its fidelity {r(t), v(t)} . 

2. For simplicity, moreover, we apply also linear-friction approximation, - replacing 
in Eqsf55] g(V) by g = g = const (with g from (1341 ). - and, besides, calculate right- 
hand side in Eqj51] under Vq = V(t Q ) = (for anyway long-time behavior of Kd(t) 
must be indifferent to V(t ) §). 

5 At that, we merely omit from solution AR of EqsJ^H] additive Vq 's contribution, exp (— gt) Vq . 




(51) 
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Then, introducing functions 
f(t) = f(t,R(t))/rn , n(t) = Vf(t,R(t))/m , C(r) = [1 - exp(-gr))/g , 
and taking t = and naturally R(0) = , we can write solution to Eqsf25] as 

AR(t) = [ C(t- t') £{t') dt' , (52) 
Jo 

while solution to Eqsj38] as infinite iteration series 

r(t) = C(t) + / C(t-t')r](t')r(t')dt' = 



C{t) + f C(t- t') r](t') C(t') dt' + (53) 
Jo 



+ / C(t- t') r](t') / C(t' - t") v {t") C(t") dt' dt" + 



Next, inserting all this into Eqj51] and imagining, - in the spirit of time- local linear 
approximation, - £(£), i](t) like white noises, one can see that in fact the only third 
term of expansion (1531) survives after averaging. It produces 

/ K D (r)dr = -(AR(t),AR(t),r(t)) w (54) 
Jo m 

dt' [ dt" C\t - t')) C(t' - t") C(t - t") C(t") x 
Jo 

J j ( I{k - t', R(h) - R{t')) I(h - t", R(t 2 ) - R{t")) ) dh dt 2 
Here new (tensor) function I(r,r) appears defined by "?7-£" cross-correlator: 

(Vf(t',R(t'))*f(t 1 ,R(t 1 ))) R /m 2 = Iitr-t^Rit^-Rit')) , 

J(r,r) = V* V* V K xx (T,r)/m 2 = (55) 



T 

rn 



x 



— - / / k*k*k sin kr cos cut S(lo, \k\) — dk 



du 

m' 2 J J 2tt 
with (. . .)r standing for conditional average under given BP's trajectory. 

Evidently (and importantly), when performing integrations over t\ and t% in Eq 
we have to make replacement 

R(ti) -R(t') =>• V(t')r - gV(t')T 2 /2 (r = ti-t'), (56) 

and similarly for R(t 2 ) — R{t") . At that, we remove from this displacements their 
parts containing f(t,R(t)) , since contribution of these parts to the average in fl54l) is 
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definitely negligible, at least under condition r v ^> r c . Accordingly, - applying under 
integral in Eqj55l approximation 

sin [k (Vt - gVr 2 /2)} « sin Wr - k (gVr 2 /2) cos Wt , 

- we can write 

J I(h - f , R(h) - R(t')) dh g T{V{t')) V{t') , (57) 

1 f IT 
T{V) = — / k 4 S"(kV,\k\)dk = —%g(V) (58) 



m 2 J m 

( u = — aJ2 — 

and analogously for second multiplier under the average (using symbolical scalar nota- 
tions, instead of tensor ones, or for simplicity taking in mind d = 1 ). 

At last, before inserting (157|) to ( 154)) let us make simplification as follows, 

(gT(V(t'))V(t')gT(V(t"))V(t")) =► g 2 T 2 ( V(t') V(t") > « (59) 

~ <7 2 T — exp ( -5- \t' - t"\) 
m 

T = y T(V)M(y)dV 

and besides notice that at t ^> ~ 1/g all functions C(. ..) in Eql5"4l. - except 
C(t' — t") , - can be replaced by constant 1/g §. Then we come to 

that is to non-decaying, infinitely long-range, diffusivity's correlation function. 

3. Let us recall that, by the K d (t) 's "microscopic" definition (H71) [2| |23| |24"| |20j 
I2T1 [26| [35] , in "macroscopic" (phenomenological) sense 

K d {t) = (D(t + r),D(t)) = (D(t + r)D{t)) - (D) 2 , 

where D{t) represents fluctuating diffusivity, and (D(t)) = D . Hence, our result ( 1601 
states that correlation function of BP's diffusivity fluctuations never decays to zero, 
as if D(t) were constant in time but randomly different from one BP's trajectory to 
another. 

Such "quasi-static" fluctuations are typical result of quantitatively rough theoreti- 
cal approaches to 1/f-type low-frequency fluctuations of diffusivity/mobility or other 
transport rates. For examples see e.g. [36j [37] . Nevertheless, nothing prevents such 



Notice also that in place of T it may be more correct to write T 2 = J T 2 (V)M(V) dV . 
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approaches from giving reasonable estimates of 1/ f-noise level, and they are rather cor- 
rect in prediction of characteristic long-range statistical scale invariance of transport 
processes pH [231 123 EH EO [22] , i.e. AR = R(t) - R(0) in the present case. 

Indeed, considering higher-order equilibrium cumulants (AR^) by means of obvi- 
ous generalization of exact relation (!5T|) . 



and again approximate expressions f l52|) -f l53|) . it is not too hard to see that (at t 3> t v ) 



with some coefficients c n which can be obtainedd from a recursive procedure. Such 
the asymptotic law detects essential non-Gaussianity of transport process and inappli- 
cability of the "law of large numbers" to it. 

At the same time, more accurate theories expectedly must lead to violation of such 
literal scale invariance and appearance of some slow-varying, logarithmic or power- 
law, factors in the Kr>{t) and higher-order cumulants fl62|) . For examples see e.g. 

[IHl IZHl [2H [201 12H [261 12H1 EOl [3H EH [37]. 

11. Variance of diffusuvity fluctuations. Discussion of the result 

1. Testing of more correct self-consistent approaches to our present problem we 
leave for future. Now, instead let us discuuss our spare but not trivial result ( 160]) . 

Its above derivation shows that formal "entry point" to the mentioned theoretical 
"narrow bottleneck" may be accounting for interplay between the exponential instabil- 
ity and friction (dissipation) both simultaneously induced by the medium (thermostat). 
Indeed, non-zero value of the integrated " rj-l; " correlation in EqESl is due to the sec- 
ond term of (151)1) reflecting BP's "braking" by the friction. Hence, in essence that 
is "instability-friction" correlation. It then naturally causes fluctuations in friction- 
related characteristics of BP's motion, first of all in diffusivity, as EqJoTl promps. 

At that, the complete " J-J" correlor in EqlHH (and its simplification in EqJ59l) is 
in fact fourth-order cross-correlator (cumulant) of the random force f(t,R(t)) and 
its gradient Vf(t,R(t)) . According to EqsJoT] [541 and finally [601 this correlation im- 
plies specific long-range fourth-order (four-point) irreducible BP's velocity correlations. 
They, in turn, in accordance with Eqj6"Tj give rise to an infinite hierarchy of higher-order 
(many-point) long-range velocity (and force) cumulants. 




(61) 



(Ai? (2n) ) « (2n- c n (AR, AR) n oc t 



(62) 
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2. Physical meaning of all such terrible picture was not once commented and ex- 
plained in our works during last thirty years (please see above references and that 
therein) . 

Our consideration once again demonstrated that diffusivity by its origin never can be 
quite certain ( a priori predictable) quantity. In absence of friction, it would have no 
definite value at all, because of constant BP's stochastic acceleration and kinetic energy 
growth (as if the medium had infinitely large temperature). But when friction takes 
place, it immediately interferes with the exponential instability (IE) and, - together 
with BP's diffusivity, - acquires fluctuations, as unboundedly diverse and unique as IE 
is in itself in (infinitely) many-particle systems. 

In principle, such consequences of IE were predicted already by N. Krylov [25] . 

3. What is for quantitative meaning of our result ( l60l) . there are many different pos- 
sibilities depending on structure of the medium's spectral function S(u, \k\) . Therefore 
let us imagine that it possesses some primitive "bell-like" shape and can be character- 
ized by three parameters only, i.e. medium's correlation time r c and length r c and, 
besides, magnitude of the force field fluctuations. The latter parameter can be replaced 
by the velocity relaxation rate g ~ 1/t v . 

Then, in view of Eqj5H]and under approximation ( 1591) . we come to estimate 

Kd r ^ ( T_ Jc_\ 2 _ ( T 

D 2 2g 2 \m r 2 J \mu 2 

where u ~ r c j r c can be interpreted as characteristic propagation velocity of the force 
field. Thus, this estimate well agrees with the reasonings expounded in Sec. 7 and 9. 

Let us also remind of, firstly, our assumption g <C l/r c , which in standard terms 
means weakness of the field fluctuations. Secondly, conventional opinion that Marco- 
vian approximation coincides with exact theory at least in the "infinitely weak inter- 
action (weak noise) limit" gr c — > . In reality, however, - as Eqj63l shows. - variance of 
relative low-frequency (quasi-static) diffusivity fluctuations, D(t)/D, is insensible to 
the gr c 's value, if the ratio r c /r c keeps constant. Hence, generally Marcovian approx- 
imation has no justification even in the weak noise limit (not speaking that anyway it 
losses any diffusivity fluctuations)! 

7 Unfortunately, conventional ideal of randomness and "stochasticity" appeals to white noise (or 
trivially related "colored noises" ) but presumes no place for something like our time-scaleless diffusiv- 
ity's (or other kinetic rates') fluctuations, even if being sowed to dynamical ground. Such an example 
is given by 39 . It shows that a biassed ideology may emasculate even mathematical physics. All the 
more so since in practice growth of scientific knowledge is accompanied by that of scientific prejudices 
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To improve our estimates of the diffusivity 1/f noise, even in case weak medium's 
noise, one should return to formally exact LE ([12!) or (1231) or may be ( 119]) and carefully 
do with their non-linearities, playing significant role at non-zero ratio r c /r c . 

12. Conclusion 

To resume, we considered Langevim equations describing random walk of particle 
in thermodynamically equilibrium fluctuating medium, and showed that the particle's 
diffusivity undergoes scaleless (1/f- type) low- frequency fluctuations whose magnitude 
can be comparable with average value of diffusivity (or even much exceed it) regardless 
of magnitude of the medium noise. 

At that we demonstrated, on one hand, usefulness of traditional "stochastic calculus" 
in the framework of dynamically based theory. On the other hand, necessity to con- 
trol stochastic way of thinking, since its seeming completeness may lead to too hasty 
and wrong conclusions (in particular it by itself automatically losses the diffusivity 
fluctuations under our interest). 

The obtained result confirms rather general "theorem on fundamental 1/f noise" 
discussed in [35] (and, under more specific conditions, in [2j|32]). I hope it will stimulate 
further investigations of transport 1/f noises in various many-particle Hamiltonian 
systems. 

(64) 
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